An acoustic analog to the dynamical Casimir effect in a Bose-Einstein condensate 
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We have realized an acoustic analog to the Dynamical Casimir effect. The density of a trapped 
Bose-Einstein condensate is modulated by changing the trap stiffness. We observe the creation of 
correlated excitations with equal and opposite momenta, and show that for a well defined modulation 
frequency, the frequency of the excitations is half that of the trap modulation frequency. 
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Although we often picture the quantum vacuum as 
containing virtual quanta whose observable effects are 
only indirect, it is a remarkable prediction of quantum 
field theory that the vacuum can generate real particles 
when boundary conditions are suddenly changed jJHi]- 
Known as the dynamical Casimir effect, a cavity with ac- 
celerating boundaries generates photon pairs. Recent ex- 
periments have demonstrated this effect in the microwave 
regime using superconducting circuits [5j [6]. Hawking 
radiation [7] is another situation characterized by spon- 
taneous pair creation and work on sonic analogs to the 
Hawking problem [8 has led to the realization that Bose- 
Einstein condensates (BEC) are attractive candidates to 
study such analog models [9HTT] , because their low tem- 
peratures promise to reveal quantum effects. Here we 
exhibit an acoustic analog to the dynamical Casimir ef- 
fect by modulating the speed of sound in a BEC. We 
show that correlated pairs of elementary excitations, 
both phonon-like and particle-like, are produced, in a 
process that formally resembles parametric down conver- 
sion ung. 

The first analyses of the dynamical Casimir effect con- 
sidered moving mirrors, but it has been suggested that 
a changing index of refraction could mimic the effect 
P~3j [14] . Our experiment is motivated by a suggestion 
in Ref. [12] that one can realize an acoustic analog to 
the dynamical Casimir effect by changing the scatter- 
ing length in an interacting Bose gas. The change in 
the interaction strength is analogous to an optical in- 
dex change: the speed of sound (or light) changes. Seen 
in a more microscopic way, the ground state of such a 
gas is the vacuum of Bogoliubov quasi-particles whose 
makeup is interaction-dependent. Changing the interac- 
tion strength projects this old vacuum onto a new state 
containing pairs of the new quasi-particles [12], which 
appear as pairwise excitations. Instead of changing the 
interaction strength, we have simply modified the confin- 
ing potential, which in turn changes the density. Sudden 
changes such as these have also been suggested as analogs 
to gravitational phenomena [15j [16] . 

We study two situations, in the first the confining po- 
tential is suddenly increased and in the second the poten- 



tial is modulated sinusoidally. The sinusoidal modulation 
of the trapping potential was studied in Ref. [T7HT9] in 
the context of the observation of Faraday waves. Our 
results on sinusoidal modulation are similar to this work 
and we have extended it to observe correlated pairs of 
Bogoliubov excitations. We produce these excitations 
in both the phonon and particle regimes, and observe 
correlations in momentum space. Parametric excitation 
of a quantum gas was also studied in optical lattices in 
which the optical lattice depth was modulated [20] [21] , 
although in that experiment, the excitation was observed 
as a broadening of a momentum distribution. 

The experimental apparatus is the same as that de- 
scribed in Refs. [22j [23] and is shown schematically in 
Fig. [pi. We start from a BEC of approximately 10 5 
metastable helium (He*) atoms evaporatively cooled in 
a vertical optical trap to a temperature of about 200 nK. 
The trapped cloud is cigar shaped with axial and radial 
frequencies of 7 Hz and 1500 Hz. In the first experiment 
we raise the trapping laser intensity by a factor of 2 with a 
time constant of 50 /is using an acousto-optic modulator 
(see inset to Fig.JTJ}). The trap frequencies thus increase 
by y/2. The compressed BEC is held for 30 ms before 
the trap laser is switched off (in less than 10 /is). The 
cloud falls onto a position sensitive, single atom detector 
which allows us to measure the atom velocities [24]. After 
compression, the gas is excited principally in the vertical 
direction: transversely we only observe a slight heating 
(about 100 nK). Figure [T|o shows a single shot distribu- 
tion of vertical atom velocities relative to the center of 
mass and integrated horizontally, while Fig. ^ shows the 
same distribution averaged over 50 shots. These distri- 
butions are more than one order of magnitude wider than 
that of an unaffected BEC. The individual shots show a 
complex structure which is not reproduced from shot to 
shot, as is seen from the washing out of the peaks upon 
averaging. 

We consider the correlations between atoms with ver- 
tical velocities v z and v' z , by constructing a normalized 
second-order correlation function, g( 2 \v z , ^)[24], aver- 
aged over the x-y plane and shown in Fig. [2}i. The plot 
exhibits two noticeable features along the v' z = v z and 
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FIG. 1: Effects of time- varying potentials (color online), a, Schematic view of the experiment. Pairs of Bogoliubov quasiparticles 
are created by varying the trap stiffness. After the flight to the detector these excitations appear as a broadening or sidebands 
on the atom cloud in the vertical (z) direction. In the following plots we convert arrival times to relative velocities and average 
over the transverse dimensions, b, Single shot velocity distribution for a cloud which was subjected to a sudden increase in the 
trap stiffness. The inset shows the time evolution of the trap stiffness, c, As in b but averaged over 50 shots, d, Single shot 
velocity distribution for a cloud which was subjected to a weak, sinusoidal modulation of the trap stiffness at 2.17 kHz. The 
inset shows the time evolution of the trap stiffness, e, As in d but averaged over 780 shots. 



v' z = — v z diagonals. The former reflects the fluctua- 
tions in the momentum distribution, as in the Hanbury 
Brown and Twiss effect [25], except that this cloud is far 
from thermal equilibrium. The v' z = — v z correlation is 
a clear signature of a correlation between quasi-particles 
of opposite velocities. A projection of this off-diagonal 
correlation is shown in Fig. [2]o. At low momentum, the 
excitations created by the perturbation are density waves 
(phonons) which in general consist of superpositions of 
several atoms traveling in opposite directions. In the con- 
ditions of our clouds, a phonon is adiabatically converted 
into a single atom of the same momentum during the re- 
lease by a process referred to as "phonon evaporation" 
[26] . Therefore in the phonon regime as well as in the 
particle regime, we interpret the back-to-back correlation 
in Fig. [2^i as the production of pairs of Bogoliubov ex- 
citations with oppositely directed momenta as predicted 
in the acoustic dynamical Casimir effect analysis [T2] . 

To further study this process, we replace the com- 
pression by a sinusoidal modulation of the laser inten- 
sity I(t) = J (l + 5cosu m t) (inset of Fig. [lji) . We 
choose 5 such that the trap frequencies are modulated 
peak to peak by about 10%. The modulation is applied 



for 25 ms before releasing the condensate. Figures[T]i and 
[]J show respectively single shot and averaged momentum 
distributions resulting from the modulation. One sees 
that the momentum distribution develops sidebands, ap- 
proximately symmetrically placed about the center. Fig- 
ure^ shows the normalized correlation function, plotted 
in the same way as in Fig. [2^i, for a modulation frequency 
co> m /27r = 2170 Hz. We again observe anti-diagonal cor- 
relations as for a sudden excitation except that the cor- 
relations now appear at a well defined velocity, which 
coincides with that of the sidebands (see Fig. ^$jp). 

We have examined sinusoidal modulation for frequen- 
cies u m /27r between 900 Hz and 5000 Hz and observed 
excitations similar to those in Fig. [3] We summarize our 
observations in Fig. [4^i in which we plot the excitation 
frequency as a function of the sideband velocity. We also 
plot the locations of the peaks in the correlation functions 
on the same graph. For modulation frequencies much 
above 2 kHz, the antidiagonal correlation functions are 
quite noisy preventing us from clearly identifying correla- 
tion peaks. This noise may have to do with the proximity 
of the parametric resonance with the transverse trap fre- 
quency (~ 3 kHz) [18]. 
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FIG. 2: Density correlations after a sudden compression 
(color online), a, Normalized correlation function g^ (v z ,v' z ) 
of clouds that have undergone a sudden compression as in 
Fig. [TJd. The signal on the diagonal results from the den- 
sity fluctuations in the cloud whereas the anti-diagonal line 
indicates the creation of correlated quasiparticles with oppo- 
site momenta. This is the signature for the analog dynamical 
Casimir effect, b, Plot of the correlation function along the 
dashed line and integrated over a region indicated by the dot- 
ted arrows, as a function of 5v z — v' z —v z . The dips on either 
side the peak may be related to the effect reported in [27] . 



A weakly interacting quantum gas obeys the well 
known Bogoliubov-de Gennes dispersion relation between 
the frequency Co>k and wavevector k: 

with a = 1 and c, the sound velocity. This relation 
describes both phonons (long wavelength excitations) 
whose dispersion is linear and free particles, whose dis- 
persion is quadratic. If our observation indeed corre- 
sponds to the creation of pairs, we expect the total exci- 
tation energy to be shared between the two excitations. 
Momentum conservation, on the other hand, requires 
that the two energies be equal, implying uj m = 2cJk- 
Therefore the relation between the modulation frequency 
and the sideband velocity should also be given by Eq. [I] 
but with a = 2 and k = mv z /h. Fitting the points in 
Fig. [4^ to with a and c as free parameters, we obtain 
a = 2.2 ± 0.3. The fitted sound velocity, 8 ± 3 mm/s, is 
consistent with the value one can calculate from the trap 
parameters and the estimated number of atoms [24] . 

We can further corroborate our interpretation of pair- 
wise excitations by a method more direct and robust than 
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FIG. 3: Density correlations after a periodic modulation 
(color online), a, Normalized correlation function (v z ,v' z ) 
measured after sinusoidal modulation of the trap frequency at 
a frequency o; m /27r = 2.17kHz, using 243 experimental shots. 
We observe a strong correlation between well defined, oppo- 
sitely directed velocities, b, Plot of the density distribution 
(blue) and of the anti-diagonal velocity correlation function, 
g {2) (v z ,v' z = -Vz) (red). 



the 2 parameter fit to the data in Fig. |4^. In Fig. [IJd, we 
compare the dispersion relation resulting from modula- 
tion with that obtained by Bragg scattering. Bragg scat- 
tering produces single excitations of quasiparticles at a 
definite energy and momentum [28 . We excited the BEC 
with two lasers in the Bragg configuration to determine 
the frequency for a given k- vector [24 . Then, under the 
same experimental conditions, using sinusoidal trap laser 
modulation, we excited the BEC at various frequencies 
and found the corresponding velocities. The lower curve 
in Fig. [4}d is a fit to the Bragg data in which we fix a = 1 
and fit the speed of sound. The upper curve is a fit to the 
trap modulation data in which we set the speed of sound 
to that found in the first fit and we allowed a to vary. 
This second fit yields a = 2.07 ±0.2. The fitted speed of 
sound for this data set (about 13 mm/s) is higher than in 
the data of Fig.[4^, because during these runs the number 
of atoms in the condensate was larger. 

An even more dramatic confirmation of our interpreta- 
tion would be the observation of sub-Poissonian intensity 
differences in the two sidebands, as was observed in the 
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FIG. 4: Dispersion relation observed by modulating the trap 
depth (color online), a, The orange squares show the side- 
band velocity determined from the density distributions. The 
green triangles are derived from the correlation functions of 
the same data. The curve is a fit to the dispersion relation 
Jl]) as described in the text. Only the solid squares were in- 
cluded in the fit: these points were all taken on the same day, 
whereas the open squares were taken under slightly different 
trap conditions, with possibly different density The error 
bars are statistical estimates based on the fits to the velocity 
distributions such as in figure ^jp. b, Comparison between 
trap modulation and Bragg scattering. The black circles are 
observations of the dispersion relation by Bragg spectroscopy. 
The orange squares are found as in (a), and clearly show that 
the corresponding frequency is about a factor of two higher 
than in the Bragg data at a given velocity The curves show 
the two fits discussed in the text. The vertical error bars on 
the Bragg data are determined by fits to the Bragg resonances. 



experiment of Ref. [5], as well as in Ref. [29] The lat- 
ter experiment modulated the center of a trapped, one 
dimensional gas producing transverse excitations which 
in turn produced twin beams. Comparing intensity dif- 
ferences in the sidebands we observe a reduction of the 
fluctuations compared to uncorrelated regions of the dis- 
tribution. However this reduction is not sub-Poissonian, 
probably because of a background under the sidebands 
(see Fig. The background is due to atoms spilling 

out of the trap before release. 

Another difference between our experiment and an 
ideal realization of the dynamical Casimir effect is that 
the temperature is not negligible. This means that the 
pair generation did not arise from the vacuum but rather 
from thermal noise. For our temperature of 200 nK, the 



thermal occupation of the mode of frequency 2 kHz is 
1.6. In the absence of a thermal background, the nor- 
malized correlation function would show an even higher 
peak. Using the perturbative approach of Ref. [12 , one 
can show that g( 2 \v Zj v' z — —v z ) is a decreasing function 
of the temperature, since thermal quasi-particles are un- 
correlated and only dilute the correlation. 

Many authors have discussed the relationship of the 
dynamical Casimir effect to Hawking radiation and Un- 
ruh radiation (see [4 for a recent review). It has also been 
noted that the two particle correlations that arise in the 
sonic Hawking problem constitute an important poten- 
tial detection technique [lOj [30] , although these authors 
examined correlations in position space. The present 
study of the dynamical Casimir effect has confirmed the 
power of correlation techniques, and shown that momen- 
tum space is a good place to look for them. We expect 
that a similar approach can be applied to Hawking radi- 
ation analogs as well as the general problem of studying 
the physics of curved spacetime by laboratory analogies. 
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Supplementary material 

Detection and data analysis. The atoms are detected 
after falling 46.5 cm to a position sensitive detector which 
allows reconstruction of the arrival time and horizontal 
position of individual atoms [25] . The mean time of flight 
is 307 ms. Given the value of the vertical Thomas- Fermi 
radius (0.5 mm), this time is long enough for the arrival 
time to reflect the vertical velocity, provided this velocity 
is well above 1.5 mm/s. 

The experimental correlation function corresponds to 
a 2D-histogram of the vertical velocities v z and v' z of 
each possible pair of atoms originating from a single con- 
densate, and averaging over all realizations. For normal- 
ization we divide this first histogram by a second 2D- 
histogram obtained with pairs of atoms originating from 
separate realizations. We observe shot to shot fluctua- 
tions in the arrival times and widths of the condensate. 
To avoid spurious correlation signals arising from fluctu- 
ations in arrival time, we recent ered each shot by aligning 
the condensate peaks before averaging and normalizing. 
The fluctuations in the width were reduced by selecting 
two different width classes and normalizing them sepa- 
rately before adding them to obtain Fig. 3. 
Speed of sound in a quasi-condensate. For our typi- 
cal atom number and trap frequencies, the atomic clouds 
are in the 1D-3D cross-over. We follow the model of 
Ref. [31] to calculate the speed of sound of atoms confined 
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in a cylindrical trap and obtain mc 2 = (jl — 
where Uj_/27r is the transverse frequency of the trap and 
fl > 1 the chemical potential in units of Huj±. We then 
apply a local density approximation for the longitudinal 
confinement on the dynamical structure factor [28] to re- 
late the speed of sound of equation (1) with the atom 
number. 

Bragg spectroscopy. In the Bragg spectroscopy mea- 
surement of Fig. 4, we use two laser beams at an an- 
gle to provide an excitation wavevector k = 2 x 
sin (0/2)k\ aser e z . The angles were approximately 6° and 
9°. The wave vector was measured by applying an in- 
tense 10 /is pulse which populated many diffraction or- 
ders. The positions of the diffraction orders permitted 
an accurate fit to find the wavevector. For Bragg diffrac- 
tion, the excitation pulse was of 5 ms duration. We then 
varied the relative frequency of the two laser beams to 
find the Bragg resonances for both positive and nega- 
tive frequency differences. The difference in the position 
of the resonances divided by 2 was used as the frequency 
co>k- The Bragg spectra exhibit a single sideband, showing 
that phonon excitations appear with a single momentum, 
as predicted by the phonon evaporation scenario. 



[2] 



[4] 



Present address: Harvard-Smithsonian Center for Astro- 
physics, Cambridge, Massachusetts 02138, USA 
t Present address: Agilent Laboratories, Santa Clara, CA 
95051, USA 

[1] G. T. Moore, J. Math. Phys. 11, 2679 (1970), ISSN 
00222488. 

S. A. Fulling and P. C. W. Davies, P. Roy. Soc. A-Math. 
Phy. 348, 393 (1976), ISSN 1364-5021. 
[3] V. V. Dodonov, Phys. Scripta 82, 038105 (2010). 

P. D. Nation, J. R. Johansson, M. P. Blencowe, and 
F. Nori, Rev. Mod. Phys. 84, 1 (2012). 
[5] C. M. Wilson, G. Johansson, A. Pourkabirian, 
M. Simoen, J. R. Johansson, T. Duty, F. Nori, and 
P. Delsing, Nature 479, 376 (2011). 
[6] P. Lahteenmaki, G. S. Paraoanu, J. Hassel, and P. J. 
Hakonen, arXiv: p. 1111.5608 (2011). 
S. W. Hawking, Nature 248, 30 (1974), ISSN 0028-0836. 
[8] W. G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 
[9] L. J. Garay, J. R. Anglin, J. I. Cirac, and P. Zoller, Phys. 
Rev. Lett. 85, 4643 (2000). 
[10] R. Balbinot, A. Fabbri, S. Fagnocchi, A. Recati, and 
I. Carusotto, Phys. Rev. A 78, 021603 (2008), URL 



jhttp : //link . aps . org/doi/10 . 1 1 03/PhysRevA . 78. ~~| 
021603l ~ 

[11] O. Lahav, A. Itah, A. Blumkin, C. Gordon, S. Rinott, 



A. Zayats, and J. Steinhauer, Phys. Rev. Lett. 105, 
240401 (2010). 

[12] I. Carusotto, R. Balbinot, A. Fabbri, and A. Recati, Eur. 

Phys. J. D 56, 391 (2010), ISSN 1434-6060. 
[13] E. Yablonovitch, Phys. Rev. Lett. 62, 1742 (1989). 
[14] F. X. Dezael and A. Lambrecht, Europhys. Lett. 89, 

14001 (2010). 

[15] P. O. Fedichev and U. R. Fischer, Phys. Rev. Lett. 91, 
240407 (2003). 

[16] P. Jain, S. Weinfurtner, M. Visser, and C. W. Gardiner, 

Phys. Rev. A 76, 033616 (2007). 
[17] P. Engels, C. Atherton, and M. A. Hoef er, Phys. Rev. 

Lett. 9 8, 095301 (2 007), URL jhttp: //link. aps. org7| 
| doi/10 . 1103/PhysRevLett . 98 . 09530l| ~ 
[18] K. Staliunas, S. Longhi, and G. J. de Valcarcel, Phys. 

Rev. A 70, 011601 (2004). 
[19] Y. Kagan and L. A. Manakova, Phys. Rev. A 76, 023601 

(2007). 

[20] C. Schori, T. Stoferle, H. Moritz, M. Kohl, and 
T. Esslinger, Phys. Rev. Lett. 93, 240402 (2004), URL 
|http : //link . aps . org/doi/10 . 1 103/PhysRevLett . 93 . | 



[21 



[22 
[23 

[24; 

[25 

[26 
[27 



240402 



C. Tozzo, M. Kramer, and F. Dalfovo, Phys. R ev. A 
72, 023613 (2005), URL |http: //link. aps.org/doi/ 10. | 
1 103/PhysRevA . 72 . 0236 13| 

G. B. Partridge, J.-C. Jaskula, M. Bonneau, D. Boiron, 
and C. I. Westbrook, Phys. Rev. A 81, 053631 (2010), 
ISSN 1050-2947. 

J.-C. Jaskula, M. Bonneau, G. B. Partridge, V. Krach- 
malnicoff, P. Deuar, K. V. Kheruntsyan, A. Aspect, 

D. Boiron, and C. I. Westbrook, Phys. Rev. Lett. 105, 
190402 (2010), ISSN 0031-9007. 

See supplementary material at 

|http: / /link.aps.org /supplementa l/XXX| 
M. Schellekens, R. Hoppeler, A. Perrin, J. V. Gomes, 
D. Boiron, A. Aspect, and C. I. Westbrook, Science 310, 
648 (2005). 

C. Tozzo and F. Dalfovo, Phys. Rev. A 69, 053606 (2004), 
ISSN 1050-2947. 

A. Imambekov, I. E. Mazets, D. S. Petrov, V. Grit- 
sev, S. Manz, S. Hofferberth, T. Schumm, E. Dern- 
ier, and J. Schmiedmayer, Phys. Rev. A 80, 
033604 (2009) , URL |http: //link. aps. org/doi/10~| 
1103/PhysRevA . 80 . 033604| 

R. Ozeri, N. Katz, J. Steinhauer, and N. Davidson, Rev. 
Mod. Phys. 77, 187 (2005), URL |http: //link. aps. org/| 



do i / 10.1 103 /RevModP hys . 77 . 187| 

R. Biicker, J. Grond, S. Manz, T. Berrada, T. Betz, 
C. Koller, U. Hohenester, T. Schumm, A. Perrin, and 
J. Schmiedmayer, Nat. Phys. 7, 608 (2011). 
P.-E. Larre, A. Recati, I. Carusotto, and N. Pavloff, Phys. 



Rev. A 85, 013621 (2012), URL http://link.aps.org/ 



doi/10. 1103/ PhysRevA . 85 . 013621]^ 

F. Gerbier, Europhys. Lett. 66, 771 (2004), ISSN 0295- 

5075. 



